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Abstract 



> 

■ A two-dimensional, fully numerical approach to the solution of four- 
(N ' 

Q^ ■ component Dirac-Fock equation using the moderately long Hermitian basis 

o 



of B-splines is applied to H, Hg" and H2 in a strong magnetic field. The geo- 



• metric parameters, including different behavior of wave- functions relativistic 

Q ■ components are analyzed. The accuracy of the solutions as a function of 

^ ■ the basis lenght is estimated. The relativistic corrections are calculated by 

transformation of the matrix equations to the equations for large relativis- 

^ • tic components. Application of the finite-element method to solution of the 

Dirac-Fock equation without supplementary assumption about exchange in 

case of the H2 excited states is discussed. 

The maximum localization of the basis functions provides applicability of 

the quadrature formulae for five-dimensional two-electron integral calculations 

within reasonable period. 

PACS number(s): 32.60.+i,33.15,03.65.Pm 
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I. INTRODUCTION 



The problem of the hydrogen atom and hydrogen molecular ion in the strong magnetic 
field has wide application in the investigations of "white dwarfs" spectral properties in 
astrophysics. Many numerical and analytical calculations have been reported in literature 
]I|-|,|-§|,|g. The moderately low temperature (5000K - 25000K) of the "white dwarfs" 
make reasonable the consideration not only of the one-electron system (H, H^), but also of 
the hydrogen molecule P-pHj], Most of the previous works treat the effect of strong 10^ 
G) magnetic fields on the electron structure in the framework of non-relativistic theory, but 
in case of the magnetic field magnitude is comparable with that of the light velocity ( 137 
in a.u.), the relativistic approach is then more consistent. First and foremost because the 
relativistic one describes the features of the particle with spin in electromagnetic field more 
correctly. In present paper we neglect the proton motion in magnetic field. Movement of 
this type may have an expressed effect on energetic structure in case of superstrong(>:^ 10^^ 
G) magnetic field ( it was estimated i.e. in p ). We aware of the fact that the relativistic 
correction in our case may be comparable to effect of the nuclear motion and this problem 
calls for further investigation. 

The one-electron system of H was studied in framework of relativistic theory in detail 
with high accuracy (see Ref. ). The molecular systems Hg" and H2 in the strong magnetic 
field were studied, as a rule, as a non-relativistic problem (see recent papers ref. [H for 
and |]9|-[ll| for H2 ). In the field- free relativistic case the maximum accuracy (better 
then 10^^° au) was obtained with use of the finite element method |[18| , ^ . The authors 
used the Lagrangian-type finite elements defined as two-dimensional fifth(and more)-order 
polynomials within triangular regions |jT3 • But for the molecular system with the exchange 



the cited authors would have to use supplementary Dirac-Fock-Slater approximation [18 



From our viewpoint the reason is that the Lagrangian-type finite elements are not enough 
localized for calculations of electron interaction matrix elements in a reasonable time (e.g. 



the fifth-order polynomial occupies 21 nodal points). As it was shown in [O] the state J2u 
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with parallel spin becomes the ground state of the hydrogen molecule in magnetic field of 
strengths B/Bq > 0.2(i?o = 2.35 x lO^G). I.e. the exchange interaction must be taken into 
account. 

The main difficulty in relativistic calculations consists in a generalized matrix problem 
enlarged by four times. Furthermore the energy matrix with an exchange loses the band 
structure and the time of calculation mainly depends on the exchange matrix construction. 
For the solution of exchange problem in a reasonable time we have used the B-splines basis 
set with minimum overlap. Taking into account the high accuracy of up-to-date relativistic 
solutions for H in the magnetic field and for in free-field case, we can estimate the 
capabilities of the finite function basis with maximum localization. The results of this 
work have shown that even with the moderate length it is possible to use this basis for 
the successful calculations of the energy relativistic corrections and geometric parameters 
of H, H2' and H2 in strong magnetic field. Following Z.Chen and S.P.Goldman |^ we have 
noted that variation method in relativistic case "may have difficulties of variational collapse, 
spurious roots, continuum dissolution in case of system with more then one electron" (in 
the one-electron system the wave functions with poor quality may have the lower-lying state 
energy, than that from the exact solution of the problem). To avoid spurious roots in case 
of H2 we used the direct QL-algorithm for finding all of the generalized eigenvalue problem 



eigenvectors |2^. Since we investigated the different behavior of one-electron wave-function 
relativistic components in the magnetic field, the calculation of H2 was performed only in 
the framework of the one-determinant approximation. 

By these means in this work we realize a two-dimensional, fully numerical approach to 
the solution of the four- component Dirac-Fock equation using the Hermitian-type basis of B- 
splines [|T^. We use the Hermitian basis set with maximum localization of each finite element 
(i.e. the overlap exists only between the nearest neighbors). In addition, the Hermitian basis 
sets provide the fast expansion and re-expansion of wave function only with use of function 
values and their derivatives in every nodal point |TB[. These features of the basis set allow 



to calculate the Dirac-Fock equation without supplementary assumption of exchange [|15 



m 



the case of the H2 states in the magnetic field. 

In this paper we have compared the energy levels of the hydrogen atom and the hydrogen 
molecular ion with the corresponding values in relativistic and non-relativistic calculations, 
have demonstrated the difference in behavior between the large and small relativistic com- 
ponents under action of magnetic field and have estimated the magnitude of magnetic field 
for transitions in the hydrogen molecule in the framework of one-determinant approach. 



II. DIRAC-FOCK FINITE-ELEMENT CALCULATION 

The finite-element method (FEM) is a well established technique for solution of several 
quantum mechanical problems. Recently, it has been demonstrated ||T7|JT8[] that this method 
can also be used for solving the many-electron Dirac problem. The finite basis functions 
presented in these publications have the following features. First, each finite function (FF) 
occupies more space than a region between the main function nodal point and the nearest 
nodal points. Second, the increasing of accuracy is provided, as a rule, by increasing the 
nodal points number for every finite function. The insufficient localization of FF impedes 
the calculation of two-electron matrix elements, therefore usually [|l7i|l8| a Slater's local 



approximation for exchange potential in diatomic molecules calculations is used. In this 
work we present a FF basis set with maximum localization for the two-center many-electron 
Dirac problem. This set we have constructed from Hermitian basis of B-splines (HBS) at 
every nodal point. The HBS analytical properties and their application to the solution of the 
one-center Fock-Dirac equations we have discussed in the above Ref. Jl^. The relativistic 



Hamiltonian operator for a diatomic molecule electron subsystem is defined as 

= i:[cc,(p,-4) + e^ft-|i-ii] + (1) 

where a and (3 are usual Dirac matrices, Zi and Z2 are the nuclear charges, ri and 12 are 
the distances to the nuclei, A = {—By, Bx, 0) is the vector-potential of the electromagnetic 
field, B is the magnetic field magnitude in z-direction. Because of axial symmetry of diatomic 
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molecules, we describe the four-component spinor of a single electron wave function $ in 
elliptic coordinates {p,T,ip): 



P 

T 



Ti + r2 

ri - r2 
D ' 



where D - is the distance between nuclei. After separation of the angular coordinate ip and 
modification of (p,t) coordinates by the introduction of the parameters for the coordinates 

"stretching" {dp and dr) 

^2 



t{T) 



P 



p- ii-d,y 

T 



(2) 



(1 + dr^ - t2 

we obtain for single electron wave function within z-projection of total momentum /i the 
following expression: 



0i(s,t) e^(''+V2)^ 
(l)2{s,t) e^(''-V2)^ 

The relativistic components 4>k are expanded in terms of the FF: 



V 



(3) 



Ms,t) = [(p^ - 1) (1 - r^)] 



2x](Ai-l/2+A)/2 



J2 C{is,iuhs,ht,k) F{is,it,hs,ht;s,t), (4) 

iB,it,hs,ht 



where {ig = l..Ng,it = l..Nt) are the nodal indices for {s,t) coordinates, {hs,ht) are the 
Hermit component indices and A=(0, 1,0,1) for k— (0,1,2,3). In our work the two-dimensional 
FF is constructed as a product of the one-dimensional elements from the HBS f{i,h;x): 



F{ig,it,hs,ht;s,t) = fiis,hs;s)f{it,ht;t) 



(5) 



The analytical properties of two first Hermit components of the B-splines f{i, h; x) are 
the following: 



f{i,0;x) 



0,\a \ > 1 
;i -a)2(l + 2a), a > 0, 
;i + a)2(l -2a), a < 0, 



(6) 



0, |a| > 1 
[1 — aY a, a > 0, 
[1 + aY a, a < 0, 

where a = {x — i d) / d, d is the distances between the nearest nodal points. 

Therefore, each FF occupies the region between the main function nodal point [ig, it) and 
the nearest nodal points 1, 1) only. We shall call this region the FF local space. We 
have calculated the H matrix elements in the crossing of the FF local spaces numerically, 
using the seven-point quadrature formula (these formulae are exact for the polynomials with 
the degree 6 and under). We have used the equi-distant net with 24 points of integration 
between the nodal points. The expansion (^) leads to a {16 NgNt) x [16 NgNt) symmetric 
eigenvalue self-consistent-field equation for the one-electron vector C 



{Hm + J - K)C = ESC , 



(7) 



where if^, J) K and S are the one-electron, Colomb's, exchange and overlap matrix op- 
erator, respectively. The number m = /i — 1/2 defines the symmetry properties of the 
wave function $. The choice of the two-dimension FF as a product of one-dimensional ones 
allows to reduce most of one-electron matrix elements to the one-dimension integrals. In 
case of the two-electron elements we used the interpolation properties of the Hermitian-type 



FF le 



Zix,y) ^ J2 C{i^,iy,h^,hy) f{i^,h^;x) f{iy,hy;y), 

'^X 



(8) 



with the nodal values 



(9) 



i.e. by the approximation of an arbitrary differentiable function the nodal values C are 
defined by the function and derivative values in the nodal points. 

Let us consider, for example, the exchange matrix element for the single-determinant 
case (in Matsuoka's notations): 

Kail =< a \ < (j)ki\g{s,t,si,ti)\al > | 0^ > = < a\Gai,ki{^,t) Icpk >, (10) 

where a = {is,it, hg, ht, k) is the generalized FF index. Function Gai,ki{s,t) is the 
result of the integration over a relatively small-sized region of the FF F{al;s,t) non-zero 
values. Moreover; in view of Eq.(|D and (P), the function Gai,ki{s,t) and its derivation can 
be calculated for the nodal points {isih) only. In the end of this section we present our 
results for Hg" ground state in the field-free case (third-order piecewise polynomials (|^)) in 
comparison with relativistic calculation on the base of fifth-order polynomials (Table |). The 
"exact" values in Table H for the relativistic orbital energies are obtained in the calculation 
of L.Yang et al. by adding the relativistic energy correction to the nonrelativistic energy 
of Madsen and Peek |^^. In Table ^ are presented the values of the ground state energy 
for H and H2 with the same parameters of FF (exept internuclear distance) as for H2 in 
Table |. The energy values for H are compared with analytical Dirac equation solution for 
Coulomb field. The calculation accuracy of H and H2 are closely allied, the accuracy for 
H2 is one order of magnitude less then for one-electron case. The relativistic corrections are 
defined by nonrelativistic limit of (|^ (a matrix equation for large relativistic components). 
The accuracy of relativistic corrections for H, H2 and H2 in the field-free case are in good 
agreement with the results of special-purpose relativistic FF calculations. The values of the 
relativistic corrections is calculated as a difference ( result of subsraction) between relativistic 
value and value in non- relativistic limit (for matrix one-electron Fock-Dirac equation (|^)) 
in the same finite functions basis. The aproach of this kind allows to compensate (in great 
part) in the relativistic corrections uniform errors, associated with the basis lenght, and to 
obtain the reasonable values with using the moderately long B-splines basis. This correction 
(as it is shown in Tables |,|| for field-free case and |T| for H in magnetic field) good fit to 
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the results, obtained with the energy high-precision calculation. 



III. RESULTS 

In the one-electron case, the equation (|^) was solved for the ground states of the hydrogen 
atom and hydrogen molecular ion. The nuclei were placed into the coordinate system focal 
points. In the coordinates (H), the Ns x A^^ equi-distant nodal net was built for the finite 
functions F from the expansion (^). The parameters dp and dr were chosen with use of the 
optimization for field-free case. The 13 x 9 for H and 10 x 11 for H2 nodal net with two 
Hermitian components of B-splines @ and four relativistic components of wave functions 
leads correspondingly to the 1872 x 1872 and 1760 x 1760 generalized matrix problem and 
to the 468 and 440 basis function per one relativistic component (See in Table | and 
Table ^ for the accuracy estimation) . This problem was solved for all of the eigenvalues by 
the method from Wilkinson and Reinsch handbook |2^ . It takes about 150 min for H or, if 
we take into account the system symmetry, 45 min for H2 of CPU time on a P-II-233 IBM 
compatible computer. The equilibrium internuclear distance was calculated with use of the 
quadratic approximation or with use of relativistic virial theorem for diatomic molecule in 
magnetic field (the specific cases of relativistic virial therem for atoms was considered in 
PJ25|1). Both of these methods lead to the near results. The results are listed in Tables |ll 



and For comparison with other results, we have used as the convention the parameter 
B/Bo, with Bo = 2.35 x lO'^G (137 a.u.). 

Our results for H presented in Table |T| are compared with detailed data from the non- 



relativistic B-splines calculations made by Jinhua Xi et al. [|T2| and with relativistic finite 
basis set calculations by Z.Chen and S.P.Goldman [^. For the field values {B/Bq from to 
10) Jinhua Xi et al. have used 35 x 16 = 560 nodal net, that is by 1.2 times more than for 
the separated relativistic component in our calculations (468). Although in our calculations 
the accuracy of energy is less then in |Q , the magnitudes of relative relativistic corrections 
are in good agreement with the results of Z.Chen and S.P.Goldman. Table shows the 
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magnetic field compression effect on the relativistic wave function too. The average radius in 
x-y-direction R^y (i.e. average radius in cyhndrical coordinat system) decreases monotonical 
both for the total density and for the small component (k=2,3 in ). But for the small 
component, one decreases to a greater extent for B/Bq values from to 20.0. The average 
length in z-direction (root-mean-square deviation of z-projection vector from nucleus) 
for the total density decreases to a greater degree then the small component. 

In the Table |I^ our results for the ground state of H2' are presented. The symmetrical 
state 1(— l/2g)^ in Herzberg's notation for double group symmetry (See |1T8| , pl| ) is considered. 

In case of antiparallel orientation of the field and z-projection of total moment (positive 
B/Bq) the relativistic correction, transversal size of the ion electronic density and size of 
small component decrease under effect of the magnetic field in the same manner as for H. 

The reverse statement is true for parallel orientation (negative B/Bq). The relativistic 
correction quadratically increases by increasing of magnetic field intensity. The similar 
behaviour can be recognized for H in Table fT| (marked with ['^]). The transversal size R^y 
of small component increases by changing B/Bq from to about -0.2, restores its original 
dimension at B/Bq ^ —1.0 and decreases at B/Bq < —1.0. 

In the Table |V|we present the "unrestricted" Dirac-Fock ( different orbitals for different 
z-projection of total momentum) results for the singlet state of II2 with the relativistic 
configuration 1(1/2^)^. The state energy has the local minimum as the function of the 
internuclear distance for all of the considered field magnitudes. The equilibrium internuclear 
distance was calculated by the quadratic approximation in the vicinity of the minimum. The 
total longitudinal and transversal sizes of the hydrogen molecule decreases, as well as for the 
hydrogen ion, but the field dependence of the small components is different for electrons with 
/i = —1/2 and fi = +1/2 (/i is z-projection of total one-electron momentum). Average 
radius of the small component in x-y-direction for the electron with fi = —1/2 decreases 
under the field compression effect monotonously. By contrast, for the electron with fi = 
+ 1/2 one increases up to the field magnitude B/Bq = 0.4 and remains greater than the 
radius in the field-free case np to B / Bq = 1.0 . The same behavior of the small component 
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is to be observed for the 1(— 1/2^) and l(+l/2g) configurations of {±B/Bq in Table |rv|) . 

The relativistic correction of total energy increases quadratically in the same manner as 
for H2 in the case of parallel orientation. 

Comparison between total energies for H2 and the sum of two energies for hydrogen 
atoms (see Table fT|) shows that aX B/Bq > 0.2 the state with two insulated atoms has the 
total energy less than that of the molecule. For better visualization of dissociation channel 
of H2 into the single atoms, we have considered the energies of the nearest configuration 
(1(— l/2g)l(— l/2u) ) with no change of the equilibrium internuclear distance. The energy 
of this configuration lie in the field-free case well above the ground state. With increasing 
of the field magnitude, this energy decreases and become equal to the 1(1/2)^ configuration 
energy in the vicinities B/Bq ^ 0.35. ..0.4. It is nessesary to stress that we have used the 
one-determinant approximation. In the non-relativistic configuration-interaction calculation 
by T.Detmer at al. |ll|] the ^J2g and ^J2u states cross over at B/Bq ^> 0.2. When taken 
into account Detmer's total energy of the ^J2g state at B/Bq = 0.2 (-1.158766a.u.) and 
total energy of two unsulated H atoms by Jinhua Xi et al. [|T2| (-1.18076313 a.u.) the cross 
point is found to be at B/Bq < 0.2. 

The energy of the 1(— l/2g)l(— 1/2^) state is free of minimum [O] as being an internu- 
clear distance function. And so, in the case of transition into 1(— l/2g)l(— 1/2„) configura- 
tion, the hydrogen molecule dissociates into single atoms. 



As it is shown in Table ^ the z-direction compression coefficient for H has given better 
fit to the relative reduction of internuclear distance in H2 and H2 than x-y compression 



one. 
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TABLES 

TABLE L Hydrogen ion ground state energies (Es), accuracy of Es, relativistic energy cor- 
rections(Rc) with increasing number of basic function per one relativistic component (Nb in our 
calculation ) or number of grid points (Np by L.Yang et al.). The values are given in au. Internu- 
clear distance for H2 is 2.0 a.u. 
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See text 
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TABLE II. H and H2 ground state energies(Es), accuracy of Es, relativistic corrections (Rc) 
with increasing number of basic function per one relativistic component (Nb). Internuclear distance 



for H2 is 


1.4 a.u. The values 


are given in 


a.u. 
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Accur.Es 


Rc 
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Accur.Es 
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24897 ^ 
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^ See |8| 
^ See Pi 
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TABLE III. Hydrogen atom ground state energies (Es), relativistic correction (Rc), relative 
correction {Rc/Enri), compression magnitude (K) for: - average radius in x-y-direction for the 
total (Rxy) density and small (Qxy) relativistic component; - dimension in z-direction for the 
total (Rz) density and small (Qz) relativistic component. In the calculations we have used 468 
of basis function per one relativistic component. The values are given in a.u.; nr - nonrelativistic 
calculations, nrl - present work nonrelativistic limit. 
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-0.831169 


nr 












1.0^ 


-0.831173226 


5.21 












2.0 


-1.0222267871 


4.031 


-4.121 


0.6072 


0.5336 


0.7716 


0.9114 


2.0 


-1.0222226658 


nrl 












2.0^ 


-1.022214 


nr 












2.0^ 


-0.022218029 


4.03 












3.0 


-1.1645485075 


3.476 


-4.048 


0.5222 


0.4634 


0.7150 


0.8738 


3.0 


-1.1645444592 


nrl 












3.0^ 


-1.164533 


nr 
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3.0^ -1.164537038 3.48 

4.0 -1.2808163924 3.117 -3.992 0.4656 0.4166 0.6752 0.8440 

4.0 -1.2808124002 nrl 

4.0'' -1.280798 nr 

8.0 -1.6194194786 2.289 -3.707 0.3463 0.3167 0.5845 0.7680 

8.0 -1.6194157712 nrl 

8.0^ -1.619384 nr 

10.0 -1.7478463122 2.015 -3.523 0.3134 0.2886 0.5574 0.7437 

10.0 -1.7478427892 nrl 

ICO"" -1.7477965 nr 

20.0 -2.2157351467 1.077 -2.388 0.2278 0.2140 0.4798 0.6714 

20.0 -2.2157327590 nrl 

20.0^^ -2.215400913 1.09 
^ Z.Chen and S.P.Goldman § 
^ Jinhua Xi et al. p| 



See comment to Table rvl 
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TABLE IV. H2 ground state electronic energies (Es), relativistic correction (Rc), total ener- 
gies(Et), equilibrium internuclear distance(Re), compression magnitude (K) of average radius in 
x-y-direction for the total (Rxy) density and small (Qxy) relativistic component. In the calcula- 
tions we have used 440 of basis function per one relativistic component. All values are given in au; 
nr - non-relativistic calculations, nrl - present work non-relativistic limit. 



7~> / 7~) a. 

B/Bq^ 


Es 


Rc 

X 10 


Et 


Re 


K(Rxy) 


K(Qxy) 


0.00'^ 




nr 


-0.60263462 


1.99719 






0.00 


-1.10334036 


nrl 


-0.60263788 


1.997194 


(0.957659) ^ 




0.00 


-1.10334774 


-7.3772 


-0.60264525 


1.997194 


1.0(0.95764) 


1.0(1.08362) 


0.02^ 




nr 


-0.61257052 


1.99699 






0.02 


-1.113327412 


nrl 


-0.61257377 


1.99697 






0.02 


-1.113334597 


-7.1850 


-0.61258096 


1.99697 


0.99985 


0.99449 


0.10^ 




nr 


-0.65103820 


1.99221 






0.10 


-1.15299661 


nrl 


-0.65104149 


1.99221 


0.99647 




0.10 


-1.15300313 


-6.5229 


-0.65104801 


1.99221 


0.99646 


0.97084 


-0.10 


-1.05300513 


-8.5243 


-0.55105001 




0.99647 


1.02410 


0.20 


-1.20187690 


-5.8942 


-0.69633157 


1.97806 


0.98651 


0.93907 


-0.20 


-1.00188103 


-10.027 


-0.49633570 




0.98652 


1.04133 


1.00^ 




nr 


-0.97321 


1.76 






1.00 


-1.54576549 


nrl 


-0.97499184 


1.752008 


0.83666 




1.00 


-1.54576963 


-4.1401 


-0.97499598 


1.752008 


0.83666 


0.72811 


-1.00 


-0.54576549 


-38.663 


0.02500815 




0.83666 


0.99836 


2.00 


-1.86877182 


-3.6984 


-1.2129584 


1.524824 


0.69635 


0.59015 


-2.00 


0.13111124 


-120.62 


0.7869246 




0.69635 


0.86803 


3.00 


-2.12263028 


-3.4913 


-1.3955526 


1.375369 


0.60793 


0.51121 


4.00 


-2.33530466 


-3.3227 


-1.5472894 


1.269011 


0.54672 


0.45855 
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-4.00 1.66425498 -443.67 2.4522702 

8.00 -2.97278069 -2.5994 -1.9998079 

10.00 -3.21982390 -2.1806 -2.17498954 

-10.0 6.77747928 -2698.9 7.82231364 
^ sign B/Bq corresponds to the field orientation 
^See Ref. [| 
^See Ref. § 

'^In the parentheses brackets the magnitudes Rxy and Qxy at B/Sq = 0.0 are shown 



0.54674 0.69928 

1.027778 0.41299 0.34670 

0.957089 0.37519 0.31551 

0.37522 0.48819 
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TABLE V. H-2 total energies(Et), electronic energies (Es), relativistic correction (Rc), equi- 
librium internuclear distance(Re) for l(l/2g)^ cofiguration of H2, compression magnitude (K) of 
average radius in x-y-direction for the total (Rxy) density and small (Qxy) relativistic component. 
In the calculations we have used 432 of basis function per one relativistic component. All values 

arc given in m.u. 





v^onng. 




xtC lie 


oneii 


xLiS 




■l^l^xyj 








xlU 
















JD / JDQ = U.U 














-i.ioooozz 


-14.0/ l.oodD 


T /'-LI /O ^ 


-u.oyo/y /y 


i.U(^i.UooD j 


1 n/' 1 ^^^^^^ 

i.ui^i.uyuoj 


1(- 


-V2ff)l(-l/2«) 


-0.7707848 


-17.17 


1/1/0 \ 

1(-V2g) 


-0.9121185 


1.0(0.9030) 


1.0(0.9969) 










-L(-J-/2uj 


n onni ctyi 


1 n/i c^nnfi\ 

i.LH^i.ouyo ) 


i.U(i.O ( /oj 








JD / Jdq U.i 














-i.izyooio 


-14. oy i.oozo 




-U.044D00O 


u.yy4Di 


n ft^i ni 

u.yoiui 










1 ('-Li /O ^ 


-U.044D0o4 


u.yy4Di 


1 nQ'?7f? 
i.Uoo(D 








B/Bq = 0.2 














-i.iloOoZo 


-io.oz i.o/iy 


1/1/0 ^ 


-U.Ooozooo 


u.yoUzu 


n noi ccn 

u.yziou 










1 / 1 1 /o ^ 
l(+l/2ff) 


-U.4ooz7i7 


u.yoUiy 


i.OooOO 








JD 1 DQ = U.OO 












1 (^ /o ^2 


1 nSQ/ISI^Q 

-i.Uoy4ooo 


1 QC 1 Q/l KQ 

-io.OO i.o4oy 




n 7/1 70 c; 1 1 

-U. <4< ZOliO 


u.youzu 


U.y44 ( 










1 / 1 1 /O "( 

l(+l/2g) 


r\ onTo^jooo 
-O.o97z6zzz 


U.98U19 


1.17371 


1(- 


-l/2,)l(-l/2„) 


-1.0493226 


-13.89 


1(-V2p) 


-1.06381193 


0.96928 


0.94771 










l(-l/2„) 


-0.33520119 


0.86675 


0.78693 








B/Bq = 0.4 












1(1/2.)^ 


-1.0768969 


-19.47 1.3373 


l(-l/2,) 


-0.7651523 


0.93671 


0.84965 










l(+l/2,) 


-0.3651650 


0.93670 


1.07798 


1(- 


-l/2,)l(-l/2„) 


-1.0806772 


-13.76 


l(-l/2,) 


-1.0827429 


0.96148 


0.88676 










l(-l/2«) 


-0.3476192 


0.84301 


0.76626 
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B/Bo = 0.8 

l(l/2g)2 -0.9370938 -33.45 1.2554 l(-l/2g) -0.8901447 0.84388 

l(+l/2g) -0.0901747 0.84387 



l(-hl/2g) 0.8459591 0.66010 



0.74181 
1.04590 



B/Bq = 1.0 

l(l/2g)2 -0.8476702 -43.69 1.21696 l(-l/2g) -0.9430515 0.80358 0.70150 

l{+l/2g) 0.0569073 0.80357 1.01730 

B/Bq = 2.0 

l(l/2g)2 -0.2894183 -126.9 1.06766 l(-l/2g) -1.1539132 0.66011 



0.57082 
0.87837 



B/Bq = 4.0 

l(l/2g)2 1.0954315 -454.2 0.898 l(-l/2g) -1.4463909 0.51404 0.44690 

l{+l/2g) 2.5531496 0.51404 0.70365 

B/Bq =10.0 

1(1/23)2 5.9491461 -2726.6 0.68682 l(-l/2g) -1.9869365 0.35056 0.31030 

l{+l/2g) 8.0103148 0.35059 0.48887 
^ In the nonrelativistic case, the configuration symbols 1(1/2^)^ and 1(— l/2g)l(— 1/2„) correspond 
Ifj^g and luglcju molecular configurations. 



See Hartree-Fock calculation by Kolos W. and Roothaan C.C.J. [26| 
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TABLE VI. The comparison of Hydrogen atom compression magnitudes K for R^y and Rz 
with relative contraction of internuclear distance D/Dq under effect of magnetic field 





K{R,) 


K{R:cy) 


D/Do 






H 


H 


Ht 




0.00 


1.0 


1.0 


1.0 


1.0 


0.10 


0.9955 


0.9914 


0.9975 


0.9970 


0.20 


0.9839 


0.9694 


0.9904 


0.9894 


1.00 


0.8632 


0.7552 


0.8772 


0.9137 


2.00 


0.7716 


0.6070 


0.7634 


0.7699 


3.00 


0.7156 


0.5222 


0.6886 




4.00 


0.6752 


0.4656 


0.6354 


0.6476 


8.00 


0.5845 


0.3463 


0.5146 




10.00 


0.5574 


0.3134 


4791 


0.4953 
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